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Z boson pole mass at two-loop order in the pure MS scheme 
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Fermi National Aecelerator Laboratory, P.O. Box 500, Batavia IL 60510 

I obtain the complex pole squared mass of the Z boson at full two-loop order 
in the Standard Model in the pure MS renormalization scheme. The input 
parameters are the running gauge couplings, the top-quark Yukawa coupling, 
the Higgs self-coupling, and the vacuum expectation value that minimizes the 
Landau gauge effective potential. The effects of non-zero Goldstone boson mass 
are resummed. Within a reasonable range of renormalization scale choices, the 
scale dependence of the computed pole mass is found to be comparable to 
the current experimental uncertainty, but the true theoretical error is likely 
somewhat larger. 
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I. INTRODUCTION 

One of the cornerstone physical observables of the Standard Model is the Z boson mass. 
The experimental value that is usually quoted is obtained using a Breit-Wigner parametriza- 
tion with a variable width, and is given in ref. [l| from a £t to LEP data as: 


= 91.1876 ±0.0021 GeV. (1.1) 

This is related j3-0] to the real part of the complex pole squared mass ~ ~ zMz 

(with P^ a constant width) according to: 


Mz = M|7P(1-P|/2M|±...) 
= 91.1535 ± 0.0021 GeV. 


( 1 . 2 ) 

(1.3) 
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In general, the complex pole squared mass is a physical observable |5-10|, independent of 
the choice of renormalization scheme and scale and the choice of gauge hxing. 

In this paper, I report a calculation, at full 2-loop order, of the complex pole squared 
mass parameters Mz and F^, using the pure MS scheme. The input parameters in this 
scheme are the running renormalized quantities 


9,9',93,yt,\v, (1.4) 

where the hrst three are the Standard Model gauge couplings, yt is the top-quark Yukawa 
coupling, A is the Higgs self-coupling, and the vacuum expectation value (VEV) v is dehned 
here to be the minimum of the full radiatively corrected effective potential in the Landau 
gauge. The normalizations used here for A and v are hxed by writing the tree-level Higgs 
potential as 


1/(4),$^) = + A(<l'^$)^ 


(1.5) 


where the canonically normalized doublet Higgs held has VEV (<h) = u/\/2, and rn? is a 
negative Higgs squared mass parameter. The minimization condition that relates v to m? 
(allowing the latter to be eliminated) is presently known at full 2-lqop order [ll| augmented 
by all 3-loop contributions at leading orders in both and yt 
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Goldstone boson mass 

effects are resummed in this relation using 13|, ; this effect is usually numerically small 

but is conceptually important, and in any case leads to simpler formulas. 

Other dehnitions of the Higgs VEV can be found in the literature. One alternative (for 


example, see refs. 15l-ll7l|) is to instead dehne the VEV as the minimum of the tree-level 
potential, so Utree = \/—vn?I\- This has the disadvantage that one must include tadpole 
diagrams explicitly. Also, one is then expanding around a point that differs from the true 
radiatively corrected vacuum, so perturbation theory converges less quickly, at least formally 
and for generic choices of the renormalization scale. Indeed, in the large yt limit, the loop 
expansion parameter is Nryf !(IQ' k'^X) rather than the usual [For more details, 

see for example refs. [13|, |l8t|, and the discussion surrounding eq. (I2.34p below.] The reason 


for the A in the denominator is that the tadpole diagrams have a Higgs propagator at zero 
momentum, which is just the reciprocal of the Higgs squared mass. 


Another alternative (see for example ref. 19||) is to dehne the VEV so that the sum of 


tadpole graphs in Feynman gauge vanishes. However, the Landau gauge effective potential is 
easier to compute to higher orders, and avoids renormalization of the gauge-hxing parameter, 
making it arguably a more convenient choice as a standard. 

The pure MS scheme is an alternative to on-shell and hybrid schemes, which have been 
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II. COMPLEX POLE MASS OF THE Z BOSON AT 2-LOOP ORDER 

To obtain the Z boson pole squared mass, one begins with the symmetric 2x2 matrix 
of neutral gauge boson transverse self-energy functions, for V, 1 /' = 7 , Z: 

nv'y'(s) = + (i 67 r 2 ) 2 ^v'y'^'®) + ••• ( 2 . 1 ) 

where s = —withp^ the external momentum, using a metric with Euclidean or 
signature. These are obtained by calculating, in the theory in d = 4—2e dimensions with bare 
parameters, the sum of 1-particle-irreducible 2-point Feynman diagrams for Ify^,, followed 
by projecting with {rj^u — p^Pu/p^)/{d — 1). The pole squared mass is then the solution of 

Spole = Zb + ^zzis^ole) + ^7^(sj^le)]V [^pole “ n 77 (sj,ie)] . ( 2 . 2 ) 


Here, 


Zb = iai + 9>i/4 


(2.3) 


is the bare, tree-level, squared mass of the Z boson. Solving eq. fl2.2p iteratively, one obtains 
to 2 -loop order: 


■^pole 




Zb + 

+in7(ZB)|Vz, 


U'^'AZb) + 


(Ibvr' 


2^12 


nW(ZB) + n«'(ZB)n«(ZB) 


(2.4) 


Instead of computing separate counterterm diagrams, the calculation described here was 
done in terms of only bare quantities qb, Qb) 93B, UtB, ^b, vb, rn?Bi then translated to 
renormalized running MS quantities g, g', g^, yti A, v at the end. Tadpole diagrams need not 
be calculated, because they automatically sum to zero, due to the dehning condition that 
the VEV is the minimum of the effective potential. Using the minimization condition for 
the Landau gauge effective potential given in ref. ISj], the parameter (and the Goldstone 


boson squared mass) are eliminated. These procedures are the same as described in refs. 18 


43| , and so most details will not be repeated here. An exception is that the 2-loop translation 
of the U{1)y gauge couplin gs f rom bare to renormalized couplings is needed, to go along 


with eqs. (2.5)-(2.24) of ref. |43| and eqs. (2.3)-(2.10) of ref. |l 8 


= h 


9 + 


" 1,1 


Ibvr^ 


+ 


(Ibvr 


2'|2 


'- 2,2 


+ 


" 2,1 


+ . . . 


(2.5) 
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where 


c? 1 = 


s' 

= -9'\ 
12^ 

1,1 

s' 

1 

CO 

'2,1 

^ V 3 

s' 

1681 , 

•2,;i 

96 ^ 


Co 1 — 


Co 9 = 


^ ,j _ n j 

72 ® 24^' 


( 2 , 6 ) 

(2.7) 

( 2 . 8 ) 


and /i is th e reg ularization scale, related to the renormalization scale Q hy = ^e'^^/dvr. As 
in refs. 1^, , the results are reduced, using the Tarasov recurrence relations ^ to a set of 

1-loop basis integrals A, B and 2-loop basis integrals T, T, U, M, following the notations 


and conventions of refs. ^,[4^. The program TSIL ^ can be used to automatically and 


efficiently evaluate these basis integrals numerically. Where possible, TSIL takes advantage 
of analytical results in terms of polylogarithms, which were given in refs. 45l-l53| . In many 


cases, analytical results for the basis integrals are not available, so TSIL ernploys Runge- 
Kutta solution of differential equations in the external momentum invariant 45| , similar to 
that suggested in ref. • 


The hnal result for the 2-loop Z boson complex pole mass can be written as: 


^pole 


— M'p' — iYzMz — Z -\- 






(Ibvr^)^ 


A 


(2),QCD 


+ A 


(2),non—QCD 
Z 


(2,9) 


In the following. 


Z 

W 

t 

h 


a-+9>Vi, 

gVji, 

vyn, 

2\v^ 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


are the tree-level MS squared masses of the Z boson, W boson, top quark, and Higgs boson, 
respectively, and the couplings of the quarks and leptons to the Z boson are: 


Qf 


= vWlZh/-Q/9'V(9' + 9't 


for / = ul, Ur, dL, dn, cr, cr, ur, where 


rpUl^ _ rpU]^ _ _ rpdl^ _ _ 

rpUu _ rj-^d^ _ rpepi _ g 

Q^l ^ ^ 2/3, 

Qdi, — QdR — -1/3, 


= 1 / 2 , 


(2.14) 


(2.15) 

(2.16) 

(2.17) 

(2.18) 
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Qe, = Qen = -1, (2.19) 

Qu, = 0 . ( 2 . 20 ) 


Also, Nc = 3, and 


Uq = Uu = Ud = riL = Tie = 3 


( 2 . 21 ) 


are the the numbers of flavors of two-component quarks and leptons of each gauge transfor¬ 
mation type, {uL,dL) and ur and du and and cr, respectively. The quantities Nc, 

uq, Hu, Ud, Ur and He are kept general in the following as a way of tagging different fermion 
contributions, although they are all equal to 3 in the Standard Model. 

The 1-loop contribution is then: 


+ al)fi{t) + Nc2au^aujif2{t) + Nc{nQ — l) 


ul 


6W 


A{W) 


+Nc{nu — l)Ou^ + d^cnQttd^ + NcUda^dn + + o-l^) + /i(0) 

/AhZ -12Z‘^ - h^' 

12fT 

AW^ AW 5Z hZ Z^ \ 




9 6hh 


1811" y 


( 2 . 22 ) 


where the fermion 1 -loop integral functions are: 




flit) 

/i(0) 

/ 2 (f) = -2tB{t,t) 


^Zi?( 0 , 0 ) + ^Z, 


(2.23) 

(2.24) 

(2.25) 


The basis integrals 5(0,0), B{h,Z), and 5(11", H"), and other integral functions 

below, are always evaluated at the external momentum invariant s = Z and renormalization 
scale Q. The bottom-quark, tan-lepton, and other fermion masses have been neglected for 
simplicity, because even at 1-loop order they make a difference of less than 1 MeV in the 
real Z pole mass. However, they can easily be restored in the 1-loop part by following the 
example of the top-quark terms in the obvious way. 

The 2-loop QCD contribution can also be written in terms of the basis integral functions 
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in a few lines: 

^^),QCD ^ ^2^_c_—^ (al^ + al^)Fi{t) + 2auj^aujiF2{t) 

+ [{nQ - l)al^ + {uu - l)al^ + nga^^ + ndaljFi{0) , (2.26) 

where: 

F,{t) = ^{Z + 

+ —-- [ {24t‘^Z - 24t^ + 20tZ^ - 8Z^)B(t, tf 

+ (32^2 - mz - 48t^)A{t)B{t, t) + (56Z - 24t + lQZ‘^lt)A{tf 
-4{t - Z){l2f - 30tZ + 7Z^)B{t, t) + (296tZ - 48^^ - lMZ‘^)A{t) 

-24t^ + 22{)t^Z - l4ltZ‘^ + 23^3 , (2.27) 

Fi(0) = -^^'M(0, 0, 0, 0, 0) - 4Z5(0, 0) - yZ, (2.28) 

F2{t) = 8t(2t — Z)M(t, t, f, t, 0) + 16t T(0, f, t) 

^ Z{4t- Z) ^ 4:tZ^)B{t, tf + {16f + 80tZ)A{t)B{t, t) 

+ {8t + UZ)A{tf + (16t3 - 200^2^ + mZfB{t, t) 

FilQt^ - imZ)A{t) + 8t^ - l4t)t^Z + 48tZ‘^ . (2.29) 


The 2-loop non-QCD contribution to the Z boson pole squared mass has the form: 


(2),non—QCD 


E T/f + Z W4‘' + Z TT’+ <=<“’• 


(2.30) 


where the list of 1-loop basis integrals is 

/« = {71(h), A{t), 7l(fT), A{Z), 5(0,0), 5(t,t), 5(h,Z), 5(fT,fT)}, (2.31) 

and the list of necessary 2-loop basis integrals is: 


{/(0,0,h), /(0,0,f), 1(0,0, IT), /(0,0,Z), /(0,h,fT), /(0,h,Z), 
J(0,t,fT), J(0,fT,Z), J(h,h,h), J(h,fT,fT), I{h,Z,Z), 

I{t,t,Z), I{W,W,Z), M(0,0,0,0,0), M(0,0,0,0,fT), M(0, 0, 0, 0, Z), 
M(0, t, 0, t, W), M(0, W, 0, W, 0), M(0, W, 0, W, t), M{h, h, Z, Z, h), 
M{h, t, Z, t, t), M{h, W, Z, W, W), M{h, Z, Z, h, Z), M{t, t, t, t, 0), 
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M{t, t, t, t, h), M{t, t, t, t, Z), M{t, W, t, W, 0), M(W, W, W, W, 0), 
M(W,W,W,W,h), M(W,W,W,W,Z), S(0,0,h), S(0,0,W), S(0,t,W), 
S(h,h,Z), S(h,t,t), S(h,W,W), S(t,t,Z), S(W,W,Z), S(Z,Z,Z), 
T(h,0,0), T(h,h,Z), T(h,t,t), T(h,W,W), T(t,0,W), T(t,h,t), 
T(t,t,Z), T(W,0,0), T(W,0,t), T(W,h,W), T(W,W,Z), T(Z,0,0), 
T(0,t,t), T(0,W,W), U(h,Z,0,0), U(h,Z,h,Z), U(h,Z,t,t), 
[/(h,Z,W,W), U(t,t,0,W), U{t,t,h,t), U{t,t,t,Z), U{W,W,0,0), 
U{W,W,0,t), UiW,W,h,W), U{W,W,W,Z), U{Z,h,h,h), 

U{Z,h,t,t), U{Z,h,W,W), U{Z,h,Z,Z)}. 


(2.32) 


The coefficients cf^ and and and are quite lengthy, so they will not be listed in 
print here. Instead, they are listed in electronic form in an ancillary hie provided with the 
arXiv source for this article, called coefficients.txt. They are ratios of polynomials of 
Z, W, t, h, and v. As usual, these coefficients are not unique, because of special identities 
that relate different basis integrals in cases where the masses are not generic. 

For each of the hve-propagator M integrals for which analytical results are not available, 
the main TSIL Runge-Kutta evaluation function TSIL_Evaluate simultaneously computes 
all of the subordinate integrals S, T, U obtained by removing one or more propagator lines. 
Therefore, only 11 calls of TSIL_Evaluate are required (in addition to the relatively fast 
evaluation of the integrals that are known in terms of polylogarithms), and in total the 
numerical computation takes well under 1 second on modern computer hardware. 

I performed a number of stringent analytical checks on the calculation, similar to those 


described for the calculations of the Higgs and W boson pole masses in 18|, l43|. First, 


’pole 


is free of poles in e. The cancellation of these poles relies on agreement between 
the divergent parts of the loop integrals performed here and the counterterm coefficients 
which can be obtained from the 2-loop scalar anomalous dimension and (3 functions from 


2„,2 


refs. [55H58[. Second, poles and logs of the Goldstone boson squared mass G = + X^v 


were checked to cancel after the resummation described in 13|, ll^. Third, I checked the 
cancellations between contributions from unphysical vector propagator components with 
poles at 0 squared mass and the corresponding Landau gauge Goldstone boson propagators. 
This ensures the absence of unphysical imaginary parts of the complex pole squared mass. 
Note that F^ = 0 in the case Nc = rii = Ue = 0. Next, I checked the absence of singularities 
in various formal limits (none of which are close to being realized in the actual parameters 
of the Standard Model), in which one or more of the following quantities vanish: Z, W, t, 
h, t — W, At — Z, AW — h, and AZ — h. This is despite the fact that many of the individual 
2 -loop coefficients do have singularities in one or more of those cases; non-trivial relations 
between basis integrals are responsible for the smooth limits of the total. Finally, I checked 
analytically that the complex pole squared mass is renormalization group scale-invariant up 














to and including all terms of 2 -loop order, using 




9 d \ d 


’pole 


= 0 , 


(2.33) 


where 7 is the Higgs anomalous dimension, and X = {(?, 5 '^ f/s, ?/*, A}. In the conventions 
used here, the derivatives of the 1 -loop basis integrals with respect to squared masses are 
listed in eqs. (A.5) and (A. 6 ) of ref. while the derivatives of the 1 -loop and 2 -loop basis 
inte gral s with respect to the renormalization scale Q can be found in eqs. (4. 7)-( 4.13) of 


ref. 451. The beta functions and scalar anomalous dimension are listed in refs. 
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In the next section, a numerical check of the Q invariance will be shown. 
In refs. 


15l . Il 6 | . a calculation of the Z boson pole mass in the pure MS scheme has 


already been given. However, unlike the present paper, they expanded around the tree-level 
dehnition of the VEV, as discussed in the Introduction above. This means that even at 
1-loop order, the results take different forms. The expression for the 1-loop pole squared 
mass contribution A^^/IOtt^ given in ea. f| 2 . 22 p above appears to differ from the result of 
eq. (B.4) of ref. 15| and eq. (B.3) of ref. [l 6 | by an amount 


—I— \-SNMit) + 3hA(h) + 12WA(W) + SIT^ + QZAiZ) + 4Z^] , (2.34) 

in the notation of the present paper. There is of course no contradiction; this merely reflects 
the difference between the tree-level contributions, which are {g'^ + g'‘^)v‘^/A in this paper and 
( 5 '^ + 5'^^)'^tree/4 In refs. [15|, ll 6 |. Note in particular the presence of 1/h oc 1/A in eq. (I2.34p : at 


loop order i, the use of the tree-level VEV results in terms proportional to 1/A^. In contrast, 
there are no A —)■ 0 singularities in the present paper. A detailed comparison would be much 
more difficult at 2-loop order, as refs. [l5|, ll 6 | also relied on doing high-order expansions in 
Z/h and Z/t and 1/4 — sin^ 6w 


III. NUMERICAL RESULTS 

Consider a benchmark set of Standard Model MS parameters dehned at the input renor¬ 
malization scale Q = = 173.34 GeV: 


g{Mt) = 0.647550, (3.1) 

g'(Mt) = 0.358521, (3.2) 

yt{Mt) = 0.93690, (3.3) 

gs{Mt) = 1.1666, (3.4) 

v{Mt) = 246.647 GeV, (3.5) 

X{Mt) = 0.12597, (3.6) 
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The gauge couplings g and g' are taken to agree with ref. 191], while yt and g^ are from 


eqs. (57) and (60) of version 4 of ref. [59 


The VEV v{Mt) and the Higgs self-coupling were 
then chosen so that Mz agrees with the central value of eq. fll.311. when computed at Q = Mz, 
and Mh agrees with the current experimental central value |60t| of = 125.09 GeV, when 
computed at Q = 160 GeV using the program SMH [61[ as described in ref. j43|. With this 
set of input parameters, one also obtains = —(92.890 GeV)^ from minimization of 

the Higgs potential using SMH at Q = Mt. In this way, the experimental measurements of 
Mz and can be used to obtain the Higgs potential parameters. The choice of Q = 160 
GeV for computing Mh was explained in ref. |43t|; at this scale the effects of top-quark 


loops in the neglected electroweak 3-loop parts should be not too large. The lower choice 
of Q = Mz for computing Mz is somewhat arbitrary. One also obtains a W boson pole 
mass of Mw = 80.329 GeV, when computed at Q = Mw, using the calculation described in 


18(1 . This translates into a Breit-Wigner mass of M^^ = 80.356 GeV, using the an alog of 


to 


eq. (II.2p above. (Somewhat coincidentally, this agrees with the value found in ref. [19 
within 1 MeV, although that calculation uses a different scheme.) 

The dependences of the computed pole mass parameters Mz and T^ on the choice of Q 
are shown in hgures 13.11 and 13.21 in various approximations. These graphs are mad e by 


62 


63 


running the input parameters g, g', yt, gs, A, and v, using their 3-loop beta functions 
from the input scale Mt to the scale Q on the horizontal axis, where is computed. In 
an idealized case that Sp^jg is computed to sufficiently high order in perturbation theory, 
Mz and Tz would be independent of Q. Therefore the Q-independence is a check on the 
calculation. I End that the calculated 2-loop value of Mz varies by only about ±2 MeV from 
its median value, over the range 70 GeV < Q < 200 GeV. Below Q = 70 GeV, the scale 
dependence is much stronger. The scale dependence is smallest for Q near 100 GeV, where 
the computed Mz has its minimum, but this does not necessarily mean that this is the best 
renormalization scale; only a higher-order calculation can reduce the theoretical uncertainty. 

With regard to the width Fz, the scale dependence of the full 2-loop result is again 
about ±2 MeV from the median value over the same range 70 GeV < Q < 200 GeV. Note 
that here, including only the QGD part of the 2-loop contribution does not actually reduce 
the scale dependence much compared to the 1-loop result. This is because most of the Q 
dependence in the width arises from the runnings of the VEV and the electroweak couplings 
of the Z boson to the fermions into which it decays, and these are independent of QGD at 
the leading (1-loop) order. The result for T^ is consistent with, and slightly lower than the 
central value of, the experimental range [l| T^ = 2.4952 ± 0.0023 GeV. Of course, there are 
much better ways to calculate Fz, because the imaginary part of the 2-loop complex pole 
mass really corresponds to only a 1-loop calculation of the Z width. (Moreover, the inclusion 
of bottom-quark mass effects, neglected above for simplicity, has a larger effect on T^ than 
on Mz, and will decrease the former by an amount of order 2 MeV due to kinematics. 
There is a significant uncertainty in estimating this reduction in the imaginary part of the 
Z complex pole mass, because of the large difference between the pole and running bottom 
quark masses.) 
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FIG. 3.1: The computed pole mass Mz of the Z boson, defined by ~ ~ as 

a function of the renormalization scale Q at which it is computed, in various approximations. 
The dotted (green) line is the tree-level result Z, the short-dashed (red) line is the 1-loop 
result, the long-dashed (blue) line is the result from the 1-loop and 2-loop QCD contribution, 
while the solid (black) line is the full 2-loop order result. The input parameters g, g',yt, gs, 
and V at the renormalization scale Q are obtained by running 3-loop renormalization group 
running, starting from eqs. (j3.ip - (l3.6p . Note that the usual Breit-Wigner mass is 0.0341 
GeV larger than the Mz shown here. 



Renormalization scale Q [GeV] 

FIG. 3.2: The computed width F^ of the Z boson, defined in terms of the complex pole 
squared mass ~ ~'i^zMz, as in Figure l3Tl The short-dashed (red) line is the 1-loop 
result, the long-dashed (blue) line is the result from the 1-loop and 2-loop QGD contribution, 
and the solid (black) line is the full 2-loop order result. 
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It is important to keep in mind that the renormalization scale dependence only provides 
a lower bound on the theoretical error. Another way of investigating the robustness of the 
calculation is to take the running top-quark squared mass t in the 1-loop part eq. fl2.22p and 
perform an expansion around an arbitrary value T that can be considered to differ from t 
by an amount that is parametrically of 1-loop order. An obvious choice is to take T to be 
the (real part of the) top-quark pole squared mass. It makes sense to do this in particular 
for the 1-Ioop contribution, because the top quark mass appears only in propagators at this 
order, not as a vertex Yukawa coupling. Expanding, one Ends: 


/i(t) = /i(T) + (t-T)[(4T-2Z)S(T,T)-4A(T)-12T + 4Z]/(4T-Z) + ..., (3.7) 
f2{t) = /2(T) + (t - T) [{2Z - 12T)5(T, T) - 4A(T) + 4T] /(4T - Z) + .... (3.8) 


I have checked that if these expansions were continued to include order {t — T)^, then the 
results for the Z pole squared mass would be nearly indistinguishable from the original 
result obtained directly from the un-expanded fi{t) and / 2 (t)- However, by instead keeping 
the expansion only at first order in {t — T) as shown, one obtains an alternative consistent 
2-loop order result for the Z pole squared mass, since t — T is to be treated as formally 
of 1-loop order. This alternative consistent 2-loop order result is numerically different, 
with the difference giving an indication of the magnitude of the error made in terminating 
perturbation theory at 2-loop order. The result of using eqs. fl3.7p and fl3.8|) compared to 
the original un-expanded fi(t) and / 2 (t) is shown in Figure [T3l We see that the alternate 
consistent 2-loop result, shown as the dashed line, has a significantly worse scale dependence, 
especially at larger Q. This suggests that the scale dependence of Mz found in the original 
calculation (the solid line) is actually accidentally small, and probably underestimates the 
theoretical error. A very similar behavior was found for the W boson mass in ref. 
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IV. OUTLOOK 


In this paper I have provided a full 2-loop calculation of the Z boson complex ^le square 
mass in the pure MS scheme, to go along with similar results for the W boson 18|] and the 
Higgs boson ^ using the same renormalization scheme and the same dehnition of the VEV. 
These calculations are an alternative to the on-shell scheme results that have been widely 
used for precision studies in the Standard Model, in which Mz instead plays the role of an 
input parameter. 

The ultimate goal should be to obtain results in which the theoretical error is very small 
compared to present and projected experimental errors. The previous section shows that 
this is certainly not obtained using just the full 2-loop calculation, as the scale dependence is 
comparable to the experimental errors, and the theoretical error is probably somewhat larger. 
There is no compelling evidence or argument that the subset of 3-loop contributions that 
are QCD and top-Yukawa enhanced will be enough to ensure the dominance of experimental 
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FIG. 3.3: A close-up of the renormalization scale dependence of the computed pole mass 
Mz- The solid line is the 2-loop result, just as in Figure IHTTl The input parameters were 
chosen so that at Q = Mz the computed pole mass agrees with the experimental central value 
Mz = 91.1535 GeV from eq. (|1.3p . The dashed line is the same, but after expanding the 
1-loop part in the MS squared mass t to linear order about the value T = (173.34 GeV)^, 
using eqs. ([321) and (13. Sp . This provides an alternate consistent 2-loop order result. The two 
approximations agree near the scale Q = 77 GeV where t = T (the top-quark running and 
pole masses coincide). Note that the Breit-Wigner mass M^'^ is 0.0341 GeV larger than the 
pole mass Mz shown here. 


errors over theoretical errors. At 2-loop order, one can see from the benchmark example 
of Figure 13.11 that the QCD contribution has a much larger scale dependence, but not a 
much larger magnitude, than the non-QCD contributions, except for smaller choices of the 
renormalization scale Q where the top-enhanced QCD corrections are big. The same thing 
was noted in the comparable results for the W boson in 1^. It is therefore reasonable to 


conclude that complete 3-loop calculations will be necessary, providing a worthy challenge 
for future work. 
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